Resonant valley filtering of massive Dirac electrons 
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Electrons in graphene, in addition to their spin, have two pseudospin degrees of freedom: sublattice 
and valley pseudospin. Valleytronics uses the valley degree of freedom as a carrier of information 
similar to the way spintronics uses electron spin. We show how a double barrier structure consisting 
of electric and vector potentials can be used to filter massive Dirac electrons based on their valley 
index. We study the resonant transmission through a finite number of barriers and we obtain the 
energy spectrum of a superlattice consisting of electric and vector potentials. When a mass term is 
included the energy bands and energy gaps at the K and K' points are different and they can be 
tuned by changing the potential. 

PACS numbers: 73.22.Pr, 73.40.Gk, 71.15.AP 



I. INTRODUCTION 

The realization of stable single-layer carbon crystals, 
called graphene, has led to an intensive investigation of 
graphene's electronic properties^. Charge carriers in a 
sheet of single-layer graphene behave like "relativistic", 
chiral massless particles with a "light speed" equal to the 
Fermi velocity and possess a gapless linear spectrum at 
the K and K' points^ ^. During the last decades there 
have been a lot of theoretical and experimental attempts 
to use the spin of the electron as a carrier of informatiorP. 
The large research interest in spintronics originates from 
the fact that it promises to be smaller, more adaptable 
and faster than today's electronic devices^. Graphene in 
addition to the spin of the electron has two more degrees 
of freedom: sublattice and valley pseudospin. Valley- 
based electronics, also known as valleytronics, uses the 
valley degree of freedom as a carrier of information sim- 
ilar to the way spintronics uses electron spin. 

Valleytronics is useful when the valley encoded infor- 
mation is preserved over long distances. Defects in the 
hexagonal lattice can cause intervalley scattering and flip 
the valley state^. This can occur around graphene edges 
or adatoms that stick to the surface of graphene. How- 
ever, in order to scatter an electron from the K valley 
to the K' valley a large transfer of momentum is needed. 
Typical disorder and Coulomb-type of scattering is un- 
able to provide this momentum and in such a case the 
valley pseudospin is a conserved quantum number in elec- 
tronic transport. This allows one to use the valley pseu- 
dospin as a carrier of information. 

It was shown recently that graphene nanoribbons with 
zigzag edges^ ^ can be used as a valley filter. Another 
promising possibility to control the valley index of elec- 
trons is by using a line of defects^^. These can be formed 
in graphene when grown on a nickel substrate or by 
using so called mass barriers that can be created by 
e.g. a proper arrangement of dopants in the graphene 
sheetP^J^. Another way to control the valley polariza- 
tion is by using local strain in graphene which induces 
an effective inhomogeneous magnetic field with opposite 



sign in the K and K' valleys^^"^. Recently, it was shown 
that a mass term can be induced by certain substrates 
such as hexagonal boron nitride (hBN) or by electron- 
electron interact ions w hich are also able to control the 
valley pseudospiiP^^^. 
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FIG. 1: (Color online) (a) Layout of the system: a ferromag- 
netic stripe on top of a graphene sheet and a hexagonal boron 
nitride (hBN) layer as a substrate, (b) Perpendicular mag- 
netic field component and corresponding vector potential at a 
distance zq from the stripe in case of parallel magnetization. 

Here we study a double Kronig-Penney (KP) model in 
graphene, consisting of a series of electric and vector po- 
tential barriers where we include a mass term as induced 
by e.g. a substrate. The vector potential can be induced 
using ferromagnetic stripes on top of a graphene layer but 
such that there is no electrical contact between graphene 
and these stripes. When one magnetizes the stripes along 
the X direction, cf. Fig. [T] by, e.g., applying an in-plane 
magnetic field, the charge carriers in the graphene layer 
feel an inhomogeneous magnetic field profile. This pro- 
file can be well approximated^^ by 2BoZoh/d{x'^ + Zq) 
on one edge of the stripe and by —2BoZoh/{x'^ + Zq) 
on the other, where zq is the distance between the two- 
dimensional electron gas (2DEG) and the stripe, and d 



and h the width and height of the stripe (see Fig. [T]). The 
resulting magnetic field profile will be modeled by two 
magnetic S functions of height 2t: B^h. Such ferromag- 
netic stripes have been deposited on top of a 2DEG in a 
semiconductor heterostructure in Ref. 22, This magnetic 
stripe can also be used as a top gate in order to create 
an electric potential barrier which will be modeled by a 
square barrier. These simplified shapes for the electric 
and magnetic field profiles will allow us to present an- 
alytical results. More realistic shapes will not influence 
the qualitative results of our work. 

The model presented in this paper is similar to a recent 
proposal^^ where a single ferromagnetic stripe on top of 
a graphene layer was shown to act as a valley filter for 
massive Dirac electrons. In our case the barrier structure 
consists of two ferromagnetic stripes acting as a resonant 
tunneling structure. We investigate both symmetric and 
antisymmetric functions of the electric and vector po- 
tentials and we show how different potential symmetries 
affect the transmission of the electrons in the K and K' 
valleys. Furthermore, a superlattice consisting of such a 
barrier structure is considered. We show that the band 
structure of the superlattice is different for K and K' 
valleys and most importantly that the height of the band 
gap is valley dependent and can be tuned by adjusting 
the potential. 

The manuscript is organized as follows. In Sec. II we 
present the model and analytical tools that were used. In 
Sec. Ill we evaluate the transmission for the combined 
electric and vector potential barrier. In Sec. IV we ex- 
pand the analysis to include multiple unit cells. In Sec. 
V we consider a superlattice of such barriers and we eval- 
uate the valley dependent band structure and band gap. 
Our concluding remarks are given in Sec. VI. 



After substituting Eq. (|3| into Eq. ([2|, we find, 

^r=E^r^ (4) 



V-\-r/jj —idx — i{ky-\-A) 



\ —idx -\- i{ky -\- A) V — Tji 

In order to find the energy spectrum for the K valley we 
use r ^ 1 and ky -^ -\-ky, while for the K^ valley we 
replace r ^ — 1 and ky -^ —ky. 



III. BARRIER STRUCTURE 

Let us first consider a single barrier consisting of both 
electric and vector potentials given by. 



v{x) = Vb [e{x) - e{x - Wb)] , 
A{x) = Ab [e{x) - e{x - Wb)] , 



(5) 
(6) 



where O is the step function, V5 and A5 are constant 
values of the electric and vector potentials and Wh is the 
width of the barrier. 

The incident wave function before the barrier is given 

by, 



i'inKX) — I Mkox 



plKQX _j_ j,p IKqX 



and the transmitted wave function after the barrier is 



(7) 



Wout[X) - y^j^^iikox+Oo) j • 



(8) 



II. MODEL 

The Dirac Hamiltonian for a massive electron in the 
presence of an external electric and vector potential is 
given hj^^^ 



H = V + rmv^cFz + vpip + eA) • a, 



(1) 



where p = —ih{dx^dy) is the momentum operator, a^ 
and a = {(Jx^cFy) are the Pauli matrices and r = ±1 is 
the valley index. The behavior of the electron is described 
by the equation. 



n^fr = E^r- 



(2) 



The equation can be simplified by using dimensionless 
units via the following substitutions: W ^ IW^ I = 1 
nm, E -^ EqE, V -^ EqV, Eq = hvp/l - 658 meV, 
A{x) -^ {Bol/p)A{x), p = eBoP/K M = rnvj,/Eo. We 
choose the Landau gauge in which A = (0, A(x)). Since 
Py commutes with the Hamiltonian of Eq. (fTl) then py 
is a good quantum number, and due to the translational 
invariance along the y direction the solutions have the 
form, 



^^rix.y) = e'^yy^ljr{x). 



(3) 



where k^ = J E'^ — fi^ — k'^ is the wave vector, with ^0 = 

arctan(A:^//co) and /o = \/E'^ — fJ? /{E + tji). For a zero 
mass term /i = the /o component reduces to sgn{E). 
Inside the barrier we have the wave function, 



i^b{x) = 



ikbX 



Afbc 



Ae 

i{kbX-\-9b) 



-ikhX 



Be 



(9) 



where /C5 = Y^^C^^^^^^P"^^7^^"^^T^y"^^^45)^ is the wave 
vector, with the angle ^5 = arctan[(/c^ + Ai))/k}j] and the 
component /^ = \/{E — VhY — l-J^'^ /{E — V5 + rpi). 

By requiring the continuity of the wave function at 
the barrier boundaries, we find the relation between the 
wave function coefficients before and after the barrier. 
This relation can be expressed as a 2 x 2 transfer matrix. 



M 



where the matrix elements are given by. 



(10) 



Mil = M*. 



22 



JkoWb 



cos{kbWb) ^ ism{kbWb) — 

\ cos ^0 cos Ob 



Mi2 = M21 = Asm6>o - smOb 

cos 6q cos ^5 



- ZTi^ COS ^0 



(11) 



(12) 



where, 



A = 



E{E - Vb) - m' 



V^2_^2^(^_y^)2_^2' 



mH 



Vi?2 - /xV(^ - ^b)" - Ai' 



(13) 



(14) 



The transmission and reflection coefficients can be found 
from the transfer matrix elements as t = l/Mn and r = 
M12/M11. 

Notice that, unlike Mn, element M12 is valley depen- 
dent. To investigate further, we can rewrite the complex 
elements of the transfer matrix to separate the amplitude 



and phase as Mn 
we find, 



miie 



M12 



mi2^ 



For Mil 



mi 



/ 2/7 TT/\ , • 2/7 ... X /sm6>osm(95- A\ 
y \ cos ^0 cos 6^5 / 

ri5 



(15) 



7 TT^ / sin 6>o sin 6>5 — A ._ tt^ nA 
^11 = k^Wb + arctan — tan^/Cfoiy^) . 

\ cos Oq cos ^6 J 

(16) 
In this case the transmission coefficient can be written as 
t = l^le*"^*. Both the transmission amplitude \t\ = l/mn 
and phase cpt = —^11 are independent of r. Next we look 

at M12, 

)^ + u'^ COS^ ^07 



sm{kbWb) / (. . n ^-^ 

^12 = 7i TTV (^sm6/o - sml95 

cos ^0 cos 6'5 



^12 = — /C0VF5 — 6>o + arctan 



Ti/ cos ^0 
A sin 6>o — sin di 



(17) 
(18) 



The reflection coefficient can be written as r = |r|e*'^'^, 
where |r| = m\2l'^\\ and (/:?r = ^\2 — ^11 • The ampli- 
tude 77112 no longer depends on r because r^ = 1. This 
means that the reflection probability i? = |rp is also in- 
dependent of the valley index. However, the phase ^\2 is 
still valley dependent. 

Electrons that are reflected at this barrier will have a 
slight phase difference depending on their valley index. 
This does not in any way affect the transmission and 
reflection probabilities for the single barrier, but, as we 
will see later, it has a significant impact on the resonance 
peaks of a double barrier structure. 



We now consider a structure that consists of two bar- 
riers as shown in Fig. [2] The left (right) barrier consists 
of an electric potential Vj, (V/^) and a vector potential 
^L (^i?)- The two barriers are separated by a distance 
Wm- 
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FIG. 2: (Color online) Double barrier structure with electric 
and vector potentials. 

The total transfer matrix M^ for the double barrier 
structure is defined as. 



Mt 



MlMmMr 



(19) 



where M^ and M^ are the transfer matrices of the left 
and right barrier, respectively, and Mm models the bar- 
rier spacing. The transfer matrix elements for the in- 
dividual barriers are the same as the previously derived 
elements in Eqs. (11) and (12) with the barrier index re- 
placed as 6 ^ L and b ^ R for the left and right barrier, 
respectively. For the middle region we define the propa- 
gation matrix Mm to connect the coefficients of the two 
barriers, 



Mm 



ikoWM 







Q ^ikoWM 



(20) 



Transmission coefficient of the double barrier structure is 
given hy t = l/M^n, so we only need to find the Mtu 
element of the total transfer matrix from Eq. (19), 



Mil = MLuMmie 



-ikWM 



MLi2MR2ie 



ikWM 



(21) 




FIG. 3: (Color online) Transmission probability for the (a) 
K and (b) K' valleys, (c) The barrier structure consists of 
a symmetric electric potential and an antisymmetric vector 
potential, (d) The difference Tr — Tk' • The solid black line 
shows the energy dispersion of a free massive electron in bulk 
graphene, E = ± y^/i^ + k^ , and the dashed lines show the 
dispersion inside the barriers, E — Vi ^ yj fi^ + (ky + AiY . 
The parameters are: Vl — Vr — 0.2, Al — —Ar — 0.1, 
/i = 0.0425, Wl^Wr^ 10, Wm = 30. 



We can substitute the complex matrix elements as Mij 
rriije^^^^ in Eq. (21) and derive the total transmissi 
probability, 

1 TlTr 



\M: 



111 



1 + RlRr + 2^/RrR^cos (t)r ' 



(22) 



where T^ (Tr) and Rl (Rr) are the transmission and 
reflection probabilities of the left (right) barrier, respec- 
tively, and (j)r is the total phase factor, 

(l)r = 2koWM - ^Lll - ^Rll + ^L12 - ^R12, (23) 

where (fLii^ ^li2 {^ru^ ^ri2) are the phase factors of 
the left (right) barrier given by Eqs. (16) and (18). Be- 
cause of the phase factor 0^- resonant transmission peaks 
will occur when. 



(l)r = (2n + l)7r, n = 0, 1,.... 

The transmission maximum is thus, 

TlTr 



(I-VWRrY' 



(24) 



(25) 



If the transmission probabilities T^ and T^ are small (as 
they are here), the reflection probabilities in the denom- 
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FIG. 4: (Color online) Transmission probability at /c^ = as 
a function of the mass term /x and the energy for the (a) K 
and (b) K' valleys, (c) The difference Tk - Tk' . The solid 
black line shows E — ii. The barrier structure is the same as 
in Fig. [3j 
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FIG. 5: (Color online) Total zero-temperature conductance 
as a function of Fermi energy for the same barrier structure 
as in Fig. |3] The inset shows the valley polarization. 



inator may be expanded as, 



Tl 



i-VSZS^.i-li-^ i-f . . ^ 



Tr\ Tl + Tr 



(26) 



CD 
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FIG. 6: (Color online) Total zero-temperature conductance as 
a function of the electric potential at Fermi energy Ep = 0.15. 
The barrier structure is the same as in Fig. [3] but with a 
variable electric potential Vl = Vr = V. The inset shows the 
valley polarization. 



and the total transmission probability is, 



[Tl+TrY 



(27) 



The minimum transmission of resonance occurs when the 
cosine function in Eq. (22) is unity, i.e. when (j)^ = 2niT. 



The single barrier transmission and reflection proba- 
bilities for the left and right barrier are completely inde- 
pendent of the valley index. As shown earlier, the phase 
ipi2 is valley dependent and as a consequence the reso- 
nant peaks of the double barrier structure will occur at 
different energies for the two valleys. 

For the valley dependent phase difference (fLi2 — ^ri2 
we find. 



^L12 

where. 



^Ri2 = -ko{WL - Wr) + arctan(7^a;), (28) 



r/i(VLAi?- VrAl), 
kofi^VLVR^QLQn 

EV,ky^{E^-^^){ky 



A,). 



(29) 
(30) 
(31) 



The valley index dependent factor is 7^- and the system 
must satisfy the condition 7^- 7^ in order to show any 
valley index dependent behavior. 



^i{VLAR - VrAl) + 0. 



(32) 



This implies that the minimum requirement is a non- 
zero mass term \i. Furthermore, both an electric and a 
vector potential are required, but this condition {VlAr — 
VrAl 7^ 0) is more flexible with regard to the individual 
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FIG. 7: (Color online) Transmission probability for the (a) 
K and (b) K' valleys, (c) The barrier structure consists of 
a symmetric vector potential and an antisymmetric electric 
potential, (d) The difference Tr — Tk' • The solid black line 
shows the energy dispersion of a free massive electron in bulk 
graphene, E = zby^/x^ + A;^, and the dashed lines show the 
dispersion inside the barriers, E — Vi ^ y^/i?~+~(S^~^h^4^. 
The parameters are: Vl — —Vr — 0.2, Al — Ar — 0.1, 
/i = 0.0425, Wl^Wr^ 10, Wm = 30. 



potential values. The system must consist of two barriers 
with different electric and/or vector potentials. In the 
case of a single barrier, i.e. Vr = and Ar = 0, condition 
( [32| cannot be fulfilled. If both a non-zero mass term 
and two barriers are present, the condition reduces to 
VlAr 7^ VrAl. Notice that the condition is still satisfied 
if one barrier has only an electric potential {Vl = V, 
Al = 0), while the other has only a vector potential 
{Vr = 0,Ar = A). 

Considering the result of the phase difference, Eq. 



(28), we can rewrite Eq. (23) as. 



(fo +r(/Pi, 



(33) 



where (po and T(pi are the valley index independent and 
dependent phase factors, respectively, 



Lpo = 2koWM - fLii - fmi 
If I = s(;n(7w)| arctan(7w)|, 
7 = KVlAr-VrAl). 



HWl 



Wr), (34) 
(35) 
(36) 



This means that the phase factor in the K valley will 
be 0x = V^o + V^i, while the K' valley will have 0^/ = 
^0 — ^1- Considering that cos (p^ is a symmetric function, 
we need to analyze the absolute value \(j)r\- The phase 
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FIG. 8: (Color online) Transmission probability at ^ = 0.28 
as a function of the mass term /i and ky for the (a) K and 
(b) K' valleys, (c) The difference Tr - Tk' • The solid black 
line shows E — ^J \i^ + k^ . The barrier structure is the same 
as in Fig. [7| 
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FIG. 9: (Color online) Total zero-temperature conductance 
as a function of Fermi energy for the same barrier structure 
as in Fig. JTl The inset shows the valley polarization. 



difference between K and K' is thus, 

^<\>^\<i>K\-\<i>K'V (37) 

If (^0 > V^i we find, 

A(/) = (/9o + ^1 - (^0 - ^i) = 2(/?i. (38) 

On the other hand if ^^ < ^\^ 

A(/) = (/9o + V^i - (^1 - ^o) = 2(/:'o. (39) 




FIG. 10: (Color online) Total zero-temperature conductance 
as a function of the electric potential at Fermi energy Ep — 
0.08. The barrier structure is the same as in Fig. [7| but with 
a variable electric potential Vl — Vr — V • The inset shows 
the valley polarization. 



Because of this, we expect to see two different resonance 
domains with different peak distributions. 

Conductance and polarization: For a given Fermi 
energy Ep^ the valley-related zero-temperature conduc- 
tance is calculated from: 



Gr = Go / Tr{EF-,ky)^ 



(40) 



where Go = e^Ly/nh and ky± = ±y^^F^~7^- The total 



conductance is G 
defined as. 



G++G_, with the valley polarization 



G+-G- 
G 



(41) 



Numerical results: The transmission probability is 
calculated numerically based on Eq. (22). We consider 
two cases. First, as shown in Fig. ^c) we take the electric 
potential of the two barriers as symmetric, while the vec- 
tor potential is taken antisymmetric. The transmission 
is shown in Figs. p[a) and (b) for K and K' , respec- 
tively. Transmission resonances are found at different 
energy values for K and K' as described by Eq. (22). 



In order to emphasize the difference we subtracted the 
electron transmission in the two valleys and plotted the 
result in Fig. [SJ^d). Notice that the difference is quite 
large and changes between -1 and 1, which means that 
for some energy values while the electron in one valley 
is completely transmitted, for the same energy we have 
complete reflection in the other valley. 

We plotted the transmission as a function of the mass 
term // and energy in Fig. |4J We can see that the two 
valleys have exactly opposite behavior with regard to the 
sign of the mass term. As the value of the mass term 







FIG. 11: (Color online) Transmission probability at /cy = 
as a function of energy and barrier spacing Wm for the (a) 
K and (b) K' valleys, (c) The barrier structure consists of 
a symmetric electric potential and an antisymmetric vector 
potential, (d) The difference Tk — Tk' ■ The parameters are: 
Vl^Vr^ 0.2, Al = -Ar = 0.1, /i = 0.0425, Wl^Wr^ 
10. 



is increased, the energy values of the resonant peaks are 
lowered for one valley and increased in the other. A no- 
table exception is the first resonant peak which behaves 
differently and vanishes into the band gap. This is a 
special case present in the domain where (/?o < ^i as 
discussed earlier. 

We evaluated the conductance and polarization using 
Eqs. ( [40l ) and ( [41] ). The results are shown in Fig. [s] 



where pronounced valley polarization is found for certain 
energies. High polarization peaks are present at lower 
values of Fermi energy, but as the energy is increased, 
the polarization becomes lower and eventually converges 
to zero. We also plotted the conductance as a function of 
the electric potential in Fig. [6J where both positive and 
negative polarization peaks are present. This shows that 
the valley filtering properties can easily be switched by 
adjusting the electric potential. 

Next, we consider the case where the electric poten- 
tial is antisymmetric, while the vector potential is sym- 
metric as shown in Fig. [tFc). The transmission in this 
case is plotted in Figs, ^a) and (b). The difference 
between the two valleys is most apparent when moving 
from ky < to ky > 0, but this is expected because of the 
ky -^ —ky substitution when switching the valleys. The 
more important result are the energy dependent trans- 
mission differences. They are a consequence of the valley 
index dependent resonances that were introduced with 
the mass term in the double barrier structure. We plot- 



ted the transmission as a function of the mass term /i 
and ky at a fixed energy in Fig. [Sj As in the previous 
case, we see that the valleys have opposite behavior with 
regard to the sign of the mass term. 

The conductance and polarization are plotted in Fig. 
|9] Compared with the previous case, the conductance is 
generally larger, but the valley polarization is much less 
pronounced. Like previously, it converges to zero for high 
energy values. Conductance is plotted as a function of the 



electric potential in Fig. 10 In this case the polarization 
is antisymmetric with respect to V = 0, which means the 
type of valley polarization can be switched by changing 
the sign of the electric potential. 

Finally, it's also possible to realize valley dependent 
transmission if the left barrier has only an electric poten- 
tial, while the right barrier has only a vector potential 
(Vl = V,Vr = 0,Al = 0,Ar = A). However, while 
interesting, there is no special advantage to this configu- 
ration. The first case, with a symmetric electric potential 
and an antisymmetric vector potential shows the largest 
energy dependent valley polarization and therefore this 
is the only configuration that will be discussed from now 
on. 

We also consider the impact of the length of the spac- 



ing between barriers Wm on the transmission. Figs. 11 
(a) and (b) show energy values of high transmission for 
the K and K' valleys and Fig. [ll[d) shows the differ- 
ence in transmission between the two valleys. The first 
transmission peak of the K' valley in Fig. pTJ^b) vanishes, 
similar to the behavior shown in Fig. Wb)"^ In this case 
the K' valley exhibits resonances from both cpo < cpi and 
(fo>^i domains, while in the K valley we always have 
(fo > (fi. 

The transmission peaks in Fig. [TTJare described by Eq. 
( |24| ). For large energy and small Wm, their positions are 
valley dependent and we can approximate them by. 




En = \hfi + U^o + ^jEl +tEA , (42) 



where, 



(2n+l)7r 



(43) 

(44) 
(45) 



On the other hand, if Wm is large, the peak positions are 
no longer valley dependent and Eq. (24) can be approx- 
imated by. 



Eo -- 


2Wt ' 


Er -- 


2^^ (Vl Vr\ 

" Wt \Al Ar) ' 


Wt -- 


= Wm + 1{Wl + Wr). 



En 




(2n + l)7r 



2W^ 



M 



(46) 



which agrees with the asymptotic behavior shown in Figs. 
11 (a) and (b). 



IV. N UNIT CELLS 

Next, we repeat the previous barrier structure TV times. 
We add a second spacing region Ws to separate unit cells 
as shown in Fig. [12] The total length of one unit cell is 
W = Wl^Wm^Wr^Ws. 

We will define a slightly different formalism than in the 
previous section. Instead of a transfer matrix that con- 
nects wave function coefficients in the different regions, 
we will define a characteristic matrix that connects the 
wave functions on the two sides of a single region^^, 

V'i(x)=M,^,(x + Wi). (47) 

The characteristic matrix for a single region is given by. 



0.24 



Mi 



1 f cos {kiWi^Oi) -ifr^sinkiWi^ 



cos 6i \ —i fi sin kiWi cos {kiWi 



(48) 



where i = L^R^M^S and the potentials in the spacing 
regions are zero: Vm = Ks = and Am = As = 0. The 
characteristic matrix for the whole structure is obtained 
by multiplying the individual matrices in the order as 
they appear in the structure M = Yl.Mi. 

For A^ unit cells the total matrix becomes M^ and 
we use a result from the theory of matrices, according 
to which the A^th power of a unimodular matrix M is 



jy _ I MiiUn-1-Un-2 M12UN-I \ ..Qx 

I M21UN-I M22UN-I - Un-2 / ' 



with X = \TrM and un is the Chebyshev polynomials 
of the second kind: 



un{x) 



sinC 



(50) 



where (" is the Bloch phase given by (" = arccos(x 
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FIG. 12: (Color online) Double barrier unit cell with electric 
and vector potentials. 
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FIG. 13: (Color online) Transmission probability and bound 
states (blue lines) for a system with 5 unit cells for the (a) 
K and (b) K' valleys. The solid black line shows the energy 
dispersion of a free massive electron in bulk graphene, E — 
± V^M^ + kl , and the dashed lines show the dispersion inside 
the barriers, E — Vi ^ yj \i^ + (ky + Ai)^. The parameters 
are: Fi, = F^ = 0.2, Al = -Ar = 0.1, /x = 0.0425, Wl = 
Wr^Wm^ 10, Ws = 5. 



The total characteristic matrix gives us the relation be- 
tween the wave functions just before and after the barrier 
structure. 



^^n(.^)=M'^^out{NW) 



(51) 



The transmission coefficient can be found by solving the 



system given by Eq. 


(51}, 


1 + r 

/o(e^^o-re-^^o) 


= te^'=«^(M^ii+MAri2/oe*^°)(52) 
= te"'°'^{MN2i + M^22/oe''°)(53) 


which leads to. 




+ — 


2g-ifcoH^cos5)p 



Mjviie- 



MnuIo + Afjv2i/o + Mjv22e»' 



(54) 

Numerical results: As an example, we consider a sys- 
tem consisting of A^ = 5 unit cells. The numerical results 
of transmission and energy states are presented in Figs. 
13 'a) and (b) for the K and K' valleys respectively. The 
valley differences are not diminished with the addition of 
more barriers. Both the transmission and bound states 
are valley dependent which suggests that a superlattice 
of this kind of system will exhibit different bands for K 
and K' . We also plot the transmission probability as a 
function of the energy E and the unit cell spacing Ws- 
As shown in Figs, p^a), (b) and (d) with an increased 
number of barriers the difference between the transmis- 
sion in K and K' is more pronounced as compared to the 
single cell case (Fig. 11). The resonances are split into 



multiple peaks due to the presence of the additional unit 
cells. When adding one extra unit cell the new electron 
wave function can be written as a symmetric or antisym- 
metric combination and each transmission peak splits to 
two. As we add A^ unit cells there are A^ eigenstates 
that arise from linear combinations of the A" states of 
the independent unit cells and as a result the maximum 
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FIG. 14: (Color online) Transmission probability at /cy = 
as a function of energy and the unit cell spacing Ws for the 
(a) K and (b) K' valleys, (c) The barrier structure consists 
of a symmetric electric potential and an antisymmetric vector 
potential, (d) The difference Tr — Tk' ■ The parameters are: 
Vl^Vr^ 0.2, Al = -Ar = 0.1, /i = 0.0425, Wl^Wm^ 
Wr = 10. 



of the transmission will split into TV peaks (e.g. in Fig. 
[Mjfor A/" = 5 we see that the resonances are split into 5 
peaks). As shown in Figs. Im] (a) and (b) for some spe- 
cific energies we have almost straight lines of full trans- 
mission. These energies correspond with the positions of 
maximum transmission shown in Figs. [tI (a) and (b) for 
and they are given by full resonant transmission 



ky 

peaks described by Eq. (|24). 



V. SUPERLATTICE 

Next, we take the same unit cell as described in the 
previous section, but consider the limit of A^ ^ oo for a 
superlattice. The spectrum of the superlattice is obtained 
from the characteristic matrix of the unit cell as, 



1 



1 



cos{k,W) = -Tr{M) = -(Mn + M22). 



(55) 



Now it is no longer meaningful to talk about transmis- 
sion. Instead, we will study the energy band structure, 
the density of state and the conductivity. 

If the spacing between unit cells Ws is much smaller 
than the width of the other regions, i.e. 1^5- <C 
Wl, Wm^ Wr^ the characteristic matrix for the cell spac- 
ing region can be replaced by an identity matrix. We will 
keep Wm proportional to Wl and Wr. The reason for 



this will be shown later. After making the substitution 
Ms -^ I and calculating the matrix for the whole struc- 
ture, the matrix elements are used in Eq. (55) to derive 



the implicit relation to obtain the energy spectrum of the 
superlattice. 



co^ikxW) = cos aL cos aM cos aR - 
jr sin aL sin aM sin aR + VLlm.R + ^m,r,l 



(56) 



where. 



r^ 



khkukR 



[AlVm^r + AmVr^l + ArVlm\ . (57) 



^H,j,k — 



rvj rC j 



£i£i 



M 



Kitvj sm ai sm aj cos a/c. 



and ai = kiWi, Vi j = Vi 



■Vj, Si 



E-Vi. 



(58) 



The spectrum of the superlattice depends on the valley 
index through 7^-. This valley dependent factor is only 
present if the superlattice has at least three potential 
regions. If we remove a region, e.g. Wm = 0, we have 
sin QfM = and 7^- is no longer a factor in the energy 
dispersion. 

A mass term is required in order to have a non-zero 
7t-. Valley dependent behavior of the superlattice also 
depends on the specific electric and vector potential val- 
ues in the three regions. This potential relation is a bit 
more complicated, but Eq. ( [57[ ) can also be represented 
in a determinant form, 



kxhkxMkxR 



Al 
Vl 



-Am 


Ar 


Vm 


Vr 


/i 


M 



(59) 



If all the values in a row or column of the determinant 
are zero then the result is zero. 

The previously derived relations are valid for any gen- 
eral case of three potentials, however we are interested 
in the specific case where the middle region corresponds 
to the space between two barriers. Therefore, in that re- 
gion, the electric and vector potentials are zero: Vm = 
and Am = 0. This will not remove 7^- completely because 
there is still a constant mass term in the second region. 
In this simpler case, the valley factor becomes. 



It 



Tfl 



kh^M^R 



{VlAr - VrAl) 



(60) 



which is similar to the valley condition Eq. ( [29| ) of the 
previously described barrier structure. We will consider 
a symmetric electric potential Vl = Vr = V with an 
asymmetric vector potential Al = —Ar = A. In that 
case the valley dependent factor becomes simply. 



-2r/i 



kLkMkR 



VA. 



(61) 
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FIG. 15: (Color online) Superlattice spectrum for the (a) K 
and (b) K' valleys. The parameters are: Vl = Vr = 0.16, 
Al = -Ar = 0.1, M = 0.0425, Wl = Wm = Wr = 10, 



FIG. 16: (Color online) Superlattice spectrum for the (a) K 
and (b) K' valleys. The parameters are: Vl = Vr = 0.29, 
Al = -Ar = 0.1, /i = 0.0425, Wl -- Wm = Wr = 10, 



It is important to mention that in all cases the val- 
ley dependent behavior requires a real external magnetic 
field. The same could not be achieved with a strain gen- 
erated field. The sign of the strain induced vector po- 
tential is different for each valley and it cancels out the 
previous valley differences. This is easy to show from 



the simplified valley factor given by Eq. (60). In the 
case of a strain induced vector potential the substitution 
Aj -^ tAj would have to be introduced. That leads to, 



r'^ji 



i^xLi^xMl^xR 



{VlAr - VrAl) ., 



(62) 



and because r^ = 1 the valley index vanishes from the 
energy dispersion relation. The same conclusion is true 
for the general three region valley factor from Eq. (57). 



Density of states (DOS): The number of k-states 



per unit energy is given by 



D{E) = -^ Yl [ dk,dkyS{E - En{k,,ky 



)). (63) 



To calculate the DOS numerically we introduce a Gaus- 
sian broadening, 



S{E-Enik^,ky)) 



1 



T^ 



exp 



{E-E^{K,ky)f 

r2 



(64) 

Conductivity: For elastic scattering the diffusive con- 
ductivity (Jij is given by 

y^ / dkxdkyTmVniVnjfn\^{l-fn\^)' (65) 



ATr'^ksT 



Here T is the temperature, v^i = dE^/dki the electron 
velocity, /nk the Fermi-Dirac function, and Tm the mo- 
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FIG. 17: (Color online) Dispersion relation (E vs kx) for ky — 
for the (a) K and (b) K' valleys with Vl^Vr^ 0.16 and 
for the (c) K and (d) K' valleys with Vl ^ Vr ^ 0.29. The 
parameters are: Al — —Ar = 0.1, /i = 0.0425, Wl — Wm — 
Wr = 10, Ws = 2. 



merit um relaxation time. For low temperatures we as- 
sume that Tm is approximately constant, evaluated at 
the Fermi level (r^ ^ tf), and replace the product 
fn\c{^-fn\^)/kBT by the delta function S{E-En{kj;, ky)). 

Numerical results: The energy bands of the su- 
perlattice are calculated numerically based on the total 
transfer matrix and the energy dispersion relation Eq. 

((55). 



As mentioned previously, we will only consider a sym- 
metric electric potential with an antisymmetric vector 
potential. The different band structures for K and K^ 
valleys with this kind of potential function are shown in 
Figs. [15] and 16 for different values of the electric poten- 
tial. In Fig. 17 we plotted the energy dispersion versus 
kx with a fixed ky = 0. In both cases the dispersion rela- 
tion is plotted for two values of the electric potential in 
order to show the large influence of the potential on the 
superlattice band structure. 

Along with the general band structure, the band gap 
depends on the valley index. Its value can be tuned by 
adjusting the electric and vector potentials, as shown in 
Fig. [Tsj At the point where either the electric or vector 
potential is zero, the band gap is identical for both valleys 
which is expected given the valley factor 7^.(61). The 
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FIG. 18: (Color online) Band gap versus the electric and vec- 
tor potentials, (a) The electric potential is tuned so that 
Vl = Vr = V, while the vector potential is taken constant 
with Al — —Ar = 0.1. (b) The vector potential is tuned so 
that Al — —Ar — A, while the vector potential is taken con- 
stant with Vl = Vr = 0.2. The rest of the SL parameters are 
the same for (a) and (b): /i = 0.0425, Wl = Wm ^ Wr -^ 10, 
Ws = 2. 



regards to the zero potential point: 

Eg{T,V,A)=Eg{-T,-V,A), 
Eg{T,V,A)=Eg{-T,V,-A), 



(66) 
(67) 



which is also obvious from the valley factor expression. 

It is possible to adjust the values so that one valley has 
a narrow band gap, while the other has a wide one, as is 
shown in Fig. 16 This case depicts the band structure at 
the point V — 0.29 from the gap function in Fig. 18 It 



is particularly convenient that the valley behavior of the 
superlattice can easily be reversed by flipping the electric 
potential from V to —F or A to —A. 

The energy spectrum and band gap of the superlattice 
also depends on the spacing between the barriers. This 
dependence is plotted in Fig. [19] For Ws — 10 there is no 
difference between K and K' . At that point the widths 
of all four regions are identical, which makes both the 
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FIG. 19: (Color online) Band gap versus the length of the 
unit cell spacing Ws for the energy spectrum of the K and 
K' valleys. The SL parameters are: Vl — Vr — 0.16, Al — 
-Ar = 0.1, M = 0.0425, Wl^Wm^Wr^ 10. 



energy gap functions for the valleys are reversed with 
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FIG. 20: (Color online) (a) Density of states (b) and (c) con- 
ductivities vs Fermi energy for K and K' valleys. The pa- 
rameters are: Vl = Vr = 0.29, Al = -Ar = 0.1, /x = 0.0425, 

Wl^Wm ^Wr = 10, Ws = 2. 



electric and vector potential functions symmetric from 
the point of view of the entire superlattice. Note, the 
potential functions are not symmetric locally in a single 
unit cell, but only globally. On the other hand, the val- 
ley difference is most pronounced at Ws = 0, because 
we have a globally symmetric electric potential with a 
globally antisymmetric vector potential. 

The numerical results for the DOS and the conductiv- 



ities a^x and dyy for V = 0.29 are plotted in Fig. 20 



These results show that the valley differences are local- 



ized to energy values around the band gap. The largest 
valley differences are present at the center of the band 
structure, while at energy values further away from the 
band gap K and K' converge to the same values, both in 
the case of the DOS and the conductivities. It should be 
noted that ayy :^ (Jxx and that cjxx goes to zero at energy 
values further away from the band gap, which means that 
the electrons have zero group velocity in the x direction. 



VI. CONCLUSIONS 

We proposed a model consisting of electric and vec- 
tor potentials for a massive Dirac electron in graphene. 
The model is essentially a series of very high and very 
narrow magnetic J- function barriers alternating in signs, 
plus electric potential barriers. The transmission through 
such a series of barriers was obtained using transfer ma- 
trix methods. 

First, we showed that although a single barrier presents 
the same transmission and reflection probabilities regard- 
less of the valley index, electrons from different valleys are 
in fact reflected with a different phase. When a second 
barrier is added, this valley phase difference will affect 
the resonant peaks that occur between the two barriers. 
We show that the transmission probability for this double 
barrier structure is very different for electrons in differ- 
ent Dirac points. This kind of structure can be formed by 
using two ferromagnetic stripes (with in-plane magneti- 
zation), which can also be used as electric gates. A mass 
term must also be present in order that this structure acts 
as a valley filter. A gap in the electronic spectrum can 
be induced either by a substrate or by electron electron 
interactions. 

We considered a few unit cells (five) and showed that 
the valley filtering behavior of the structure is robust. 
This enabled us to extend the calculation to a super- 
lattice for which the unit cell corresponds to two stripe 
ferromagnetic gates. We show that depending on the con- 
figuration of the magnetization and electric potential of 
the gates, different band structures appear for electrons 
in the two Dirac points. The generated band gap is also 
valley dependent and can be easily tuned by changing the 
electric potential or the magnetization. Even more, the 
behavior of electrons in the two valleys can be switched 
by flipping the sign of either potential (V -^ —V or 
A^-A). 

From an experimental point of view the valley filters 
and polarizers could be easily realized as superlattices (or 
even finite number of unit cells) made of ferromagnetic 
stripes all magnetized in-plane and in the same direc- 
tion but with alternating electric potentials. Alterna- 
tively, although less feasible experimentally, one can also 
obtain valley filtering behavior for alternating magneti- 
zation but with the same electric potentials. An essen- 
tial ingredient is the presence of a finite gap in the elec- 
tronic spectrum, which is nowadays routinely achieved in 
graphene deposited on BN substrates. 
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